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Abstract 

Classical matching theory can be defined in terms of matrices with 
nonnegative entries. The notion of Positive operator , central in 
Quantum Theory , is a natural generalization of matrices with non- 
negative entries. Based on this point of view , we introduce a defi- 
nition of perfect Quantum (operator) matching . We show that the 
new notion inherits many "classical" properties , but not all of them 
. This new notion goes somewhere beyound matroids . For separa- 
ble bipartite quantum states this new notion coinsides with the full 
rank property of the intersection of two corresponding geometric 
matroids . In the classical situation , permanents are naturally asso- 
ciated with perfects matchings. We introduce an analog of perma- 
nents for positive operators, called Quantum Permanent and show 
how this generalization of the permanent is related to the Quantum 
Entanglement. Besides many other things , Quantum Permanents 
provide new rational inequalities necessary for the separability of 
bipartite quantum states . Using Quantum Permanents , we give 
deterministic poly-time algorithm to solve Hidden Matroids Inter- 
section Problem and indicate some "classical" complexity difficul- 
ties associated with the Quantum Entanglement. Finally , we prove 
that the weak membership problem for the convex set of separable 
bipartite density matrices is NP-HARD. 

1 Introduction and Main Definitions 

The (classical) Matching Theory is an important , well studied but 
still very active part of the Graph Theory (Combinatorics) . The 
Quantum Entanglement is one of the central topics in Quantum In- 
formation Theory . We quote from [ pip : "An understanding of 
entanglement seems to be at the heart of theories of quantum com- 
putations and quantum cryptography , as it has been at the heart 
of quantum mechanics itself . " We will introduce in this paper 
a Quantum generalization of the Matching Theory and will show 
that this generalization gives new and surprising insights on the na- 
ture of the Quantum Entanglement . Of course , there already exist 
several "bipartite" generalizations of (classical) bipartite matching 
theory . The most relevant to our paper is the Theory of Matroids 
, namely its part analyzing properties of intersections of two geo- 
metric matroids . 

Definition 1.1 Intersection of two geometric matroids MI(X, Y) = 



{(xi, yi), 1 < i < K} is a finite family of distinct 2-tuples of non- 
zero Af-dimensional complex vectors , i.e. Xi, y, £ C . 
The rank of MI(X, Y) is the largest integer m such that there ex- 
ist 1 < ii < ... < im < K with both sets {xt 1 , Xi m } and 
{yi 1 , ...jyi^} being linearly independent. If Rank(MI(X,Y)) 
is equal to N then MI(X, Y) is called matching . The matroidal 
permanent MPrx Y) is defined as follows : 



MP, 



(X,Y) 



<i 2 <-.<iN<K 



det (Ei< fc <iv x ^ x \ k ) det (Ei< fc <iv vikVijn 



Remark 1.2: Let us denote linear space (over complex numbers 
) of N x N complex matrices as M(N) . It is clear from this 
definition that MI(X, Y) is matching iff MP( X ,y) > 0. More- 
over , MI(X, Y) is matching iff the linear subspace Lin(X, Y) C 
M(N) generated by the matrices {xiy\, 1 < i < K} contains a 
nonsingular matrix and , in general , Rank(MI(X,Y)) is equal 
to the maximal matrix rank achieved in Lin(X, Y) . The follow- 
ing equality generalizes Barvinok's (|[lo|] ) unbiased estimator for 
mixed discriminants : 



MP (x , y) = E(\ det( ti*iVi)\ 2 ) 

Ki<K 



(2) 



where 1 < i < K} are zero mean independent (or even 
2,/V-wise independent ) complex valued random variables such that 
E(\£,i\ 2 = 1, I < i < K . It is not clear whether the analysis from 



j can be applied to A/P ( 



(X,Y) 



Example 1.3: Suppose that Xi £ {ei, ejv}, 1 < i < K , where 
{ei, bat} is a standard basis in C . Define the following posi- 
tive semidefinite N x N matrices : 

(e ilBi )e(x,r) 

Then it is easy to see that in this case matroidal permanent coinsides 
with the mixed discriminant , i.e. MP(x,Y) = M(Qi, •• ■ ,Qn) 
where the mixed discriminant defined as follows : 



M(Qi,...Qn) = 



0" 



■ detfziQi + .... + x N Q N ). (3) 



dxi...dx n 

We will also use the following equivalent definition : 



M(Qi,...Q N )= (-l) Sifl " (CTr) n Qi ^W' T W)' (4) 



where S„ is the symmetric group, i.e. the group of all permutations 
of the set {1, 2, • • • , TV}. If matrices Qi,l < i < N are diagonal 
then their mixed discriminant is equal to the corresponding perma- 
nent ([§]). | 

Let us pose , before moving to Quantum generalizations , the fol- 
lowing "classical" desision problem . We will call it Hidden Ma- 
troids Intersection Problem (HMIP ) : 

Problem 1.4: Given linear subspace L C M(N) and a promise 
that L has a ( hidden ) basis consisting of rank one matrices. Is 
there exists poly-time deterministic algorithm to decide whether L 
contains a nonsingular matrix ? Or more generally , to compute 
maximum matrix rank achieved in L ? | 

Below in the paper we will assume that linear subspace L C M(N) 
in (HMIP ) is given as a some rational basis in it. If this basis con- 
sists of rank one matrices then there is nothing "hidden" and one 
can just apply standard poly-time deterministic algorithm comput- 
ing rank of intersection of two matroids. A natural (trivial) way 
to attack (HMIP ) would be to exract a ( hidden ) basis consisting 
of rank one matrices. We are not aware about the complexity of 
this extraction . The following example shows that there exist lin- 
ear subspaces L C M(N) having a rational real basis and a "rank 
one" basis but without rational "rank one" basis : 
Consider the following 2x2 matrix 



.4 



-2 
1 



and define linear subspace IR C M(2) generated by A and the 

identity I . 

It is easy to see that 

Rank{aA + bl) < 1 iff a 2 + 2b 2 = 0. 

Therefore there are no rank one rational (complex) matrices in IR 
. From the other hand rank one matrices 

C = V2I + iB,D = V2I - iB 

form a basis in IR . 

One of the main results of our paper is a positive answer to the 
nonsingularity part of (HMIP ) . Moreover our algorithm is rather 
simple and does not require to work with algebraic numbers . 
And , of course , we are aware about randomized poly-time algo- 
rithms , based on Scwartz's lemma , to solve this part of (HMIP 
) . But for general linear subspaces , i.e. without extra promise , 
poly-time deterministic algorithms are not known and the problem 
is believed to be "HARD" . To move to Quantum generalization , 
we need to recall several , standard in Quantum Information litera- 
ture , notions . 

1.1 Positive and completely positive operators ; bipartite 
density matrices and Quantum Entanglement 

Definition 1.5 A positive semi definite matrix pa,b '■ C N ®C N — > 
C N C N is called bipartite unnormalized density matrix 
(BUDM ) , if tr(pA.B) = 1 then this Pa,b is called bipartite den- 
sity matrix . 

It is convinient to represent bipartite pa,b = p(ii, *2, 31,32) as the 
following block matrix : 



Pa,b = 



( Ax,i 


Al t 2 ■ 


■ Al,N 


A 2 A 


A 2> 2 ■ 


A?,N 




An, 2 


■ An,n 



(5) 



where A iujl =: {p(ii,ii,ji,ja) ■ 1 < i2,32 < N},1 < 
< TV. 

A (BUDM ) p called separable if 



P = P(X,Y) =■ ^2 XiX i ® ViV *> 
l<i<K 



(6) 



and entangled otherwise . 

If vectors Xi , yi ; 1 < i < if in (6) are real then p is called real 
separable . 

Quantum marginals defined as pa = ^2 1<i<N A%,i and 
(p B (i,j)=tr(Aijy,l<i,j < N) . 

We will call (BUDM ) p weakly separable if there exists a sep- 
arable p[ x Y \ with the same Image as p : Im(p) = Im(p'^ x Y ))- 
A linear operator T : M(N) -> M(N) called positive if T(X) y 
for all X t_ , and strictly positive ifT(X) y atr(X)I for all 
X y and some a > 0. A positive operator T is called completely 
positive if 



T{X) 



E 

Ki<N 2 



AiXA\;Ai,X G M(N) 



(7) 



Choi's representation of linear operator T : M(N) — * M(N) is a 
block matrix CH(T)ij =: T(eiet). Dual to T respect to the inner 

product < X, Y >= tr(XY' 1 ) is denoted as T*. Very usefull and 
easy Choi's result states that T is completely positive iff CH(T) 
is (BUDM ) . Using this natural (linear) correspondence between 
completely positive operators and (BUDM ) , we will freely "trans- 
fer" properties of (BUDM ) to completely positive operators . For 
example , a linear operator T is called separable iff CH(T) is sep- 
arable , i.e. 

r(Z)=T (JflY )(Z)= **Vi Z Vi*i ( § ) 

l<i<K 

Notice that CH(T {X>Y) ) = p iX)Y) and T ( * x Y) = T {Y ,x) ■ I 

Remark 1.6: In light of definition (1.5) , we will represent linear 
subspaces L C M(N) ^ C N ® C N in (HMIP ) as images of 
weakly separable (BUDM ) p . And as the complexity measure we 
will use the number of bits of (rational) entries of p. | 

The next definition introduces the quantum permanent QP(p) , the 
main tool to solve (HMIP ) . Though it was not our original inten- 
tion , it happens that QP(pix,Y)) = MP(x,y) ■ 

Definition 1.7: We define quantum permanent, QP(p) , by the 
following equivalent formulas : 



QP(p) 



E 



(-ir 3n{a) M(A 1M1) ,...,A NMN) ); (9) 



QP(P)= ]T (-ir 9 " (T1T2T3) n^,n(i),r2(i),r 3 (i)); 

Tl,T2,T 3 £Sx 1=1 

(10) 



= m E 

N 

J|r/ l0 (T 1 (i),T 2 ( l ),r 3 ( l ),r4(j)). (11) 



I 



Remark 1.8: The representation (6) is not unique , it follows di- 
rectly from the Caratheodory Theorem that one always can choose 
K < N 4 in (6) . Thus , the set of separable (BUDM ) , de- 
noted by Sep(N, N) , is a convex closed set . As it is known that 
Sep(N, N) has non-empty interiour , it follows from straightfor- 
ward dimensions counting that for the "most" separable (BUDM ) 
at least K > ^frr- 1 

In the next proposition we summarize the properties of the quantum 
permanents we will need later in the paper . 

Proposition 1.9: 

1. 

QP(P(X,Y)) = MP {X , Y ) (12) 

2. 

QP(p)=<p® N Z,Z>, (13) 

where p® N stands for a tensor product of N copies of p , 
< . , . > is a standard inner product and 

if3k } = Tfc(i)(l < i < N);r k G S N (k = 1,2) and zero 
otherwise . 

( The equality (13) implies that if pi X P2 tl then 
QP{pi)>QP{p 2 )>0.) 

3. 

QP{(A! ® A 2 )p(A 3 ® A 4 ) = det(A 1 A 2 A 3 A4)QP(p) 

(14) 

4. 

QP{pa,b) = QP(pb,a) (15) 

Example 1.10: Let us present a few cases when Quantum Perma- 
nents can be computed "exactly ". They will also illustrate how 
universal is this new notion . 

1. Let pa,b be a product state , i.e. pa.b — C C8> D . Then 

QP{C® D) = Det(C)Det(D) . 

2. Let Pa,b be a pure state , i.e. there exists a matrix (R = 

R(i,j):l<i,3<N) 

such that pA,B{ii,i2,h,h) = R{h, i2)R(ji, ja) ■ 
In this case QP{pa,b) = N\\Det(R)\ 2 . 

3. Define blocks of pa.b as Aij = R(i,j)eie\ . 
Then QP{pa,b) = Per(R) . 

I 

The next definition introduces Quantum Perfect Matching. 

Definition 1.11: Let us consider a positive (linear) operator T : 
M(N) -> M(N) , a map G : C — > C N , and the following 
three conditions : 

1. G(x) G Im(T{xx^). 

2. If {xi, xn} is abasis in C N then {G(xi), G(xn)} is 
also a basis, i.e. the map G preserves linear independence. 

3. If {xi, rcjv} is an orthogonal basis in C N then 
{G(xi), G(xn)} is a basis . 

. We say that map G is Quantum Perfect Matching for T if it sat- 
isfies conditions ( 1 ,2) above ; say map G is Quantum Semi-Perfect 
Matching for T if it satisfies conditions (1,3) above . 
I 



In the rest of the paper we will address the following topics : 

1 . Characterization of Quantum Perfect Matchings in spirits of 
Hall's theorem . 

2. Topological and algebraic properties of Quantum Perfect Match- 
ings , i.e. properties of maps G in Definition (1.11). 

3. Compelexity of checking whether given positive operator is 
matching . 

4. Quantum (or Operator ) generalizations of Sinkhorn's itera- 
tions (in the spirit of [Q , |§§] , Q ). 

5. van der Waerden Conjecture for Quantum Permanents. 

6. Connections between topics above and the Quantum Entan- 
lement . 

7. Complexity to check the separability . 

2 Necessary and sufficient conditions for Quantum Perfect 
Matchings 

Definition 2.1: A positive linear operator T : M(N) — » M(N) 
called rank non-decreasing iff 

Rank(T(X)) > Rank(X) ifX y 0; (16) 

and called indecomposable iff 

Rank{T(X)) > Rank(X) ifX V andl < Rank(X) < N. 

(17) 

A positive linear operator T : M(N) — > M(N) called doubly 
stochastic iff T{I) = / and T*(I) = I ; called e - doubly stochas- 
tic iff DS{T) =: tr{(T{I) - I) 2 ) + tr({T*{I) - I) 2 ) < e 2 . 
I 

The next conjectures generalize Hall's theorem to Quantum Perfect 
Matchings . 

Conjecture 2.2: Assuming that the Axiom of Choice and the Con- 
tinium Hypothesis hold, a positive linear operator T has Quantum 
Perfect Matching iff it is rank non-decreasing . | 

Conjecture 2.3: Assuming that the Axiom of Choice and the Con- 
tinium Hypothesis hold, a positive linear operator T has Quantum 
Semi-Perfect Matching iff it is rank non-decreasing . | 

Remark 2.4: We realize that the presence of the Axiom of Choice 
and the Continium Hypothesis in linear finite dimensional result 
might look a bit weird . But we will illustrate below in this section 
that for some completely positive entangled operators correspond- 
ing Quantum semi-perfect matching maps G are necessary quite 
complicated , for instance necessary discontinues . Moreover Con- 
jecture 1 is plain wrong , even for doubly stochastic indecompos- 
able completely positive operators . In separable and even weakly 
separable cases one does not need "exotic axioms" and one can re- 
alize Quantum perfect matching map it it exists as a linear nonsin- 
gular transformation through a rather simple use of Edmonds-Rado 
theorem . | 

The next Proposition(2.5) is a slight generalization of the corre- 
sponding result in | |24| ] . 

Proposition 2.S Doubly stochastic operators are rank non-decreasing 
. If either T(I) = I or T*{I) = I and DS(T) < N' 1 then T 
is rank non- decreasing . IfDS(T) < (2N + 1) _1 then T is rank 
non-decreasing . 



Example 2.6: Consider the following completely positive doubly 
stochastic operator Sk 3 : M(3) -> M(3) : 

Sk :i (X) = iA (1 , 2 )A-Aj li2) + J 4 (1 ,3)XA (li3) t+yl(2,3)^(2,3)t 

(18) 

Here {A( it j), 1 < i < j < 3} is a standard basis in a linear 
subspace of M(3) consisting of all skew- symmetric matrices , i.e. 

=: Cie] — e^ej and {e;,l < i < 3} is a standard or- 
thonormal basis in C 3 . It is easy to see that for a real normed 
3-dimensional column vector x the image ImSks(xx') is equal to 
the real orthogonal complement of x , i.e. to the linear 2-dimensional 
subspace x of R 3 consisting of all real vectors orthogonal to x 
. Suppose that G is Quantum semi-perfect matching map , then 
G(x) G x ± and , at least , G(x) is nonzero for nonzero vectors 
x. By the well known topological result , impossibility to comb the 
unit sphere in R 3 , none of Quantum semi- perfect matchings for 
Skz is continuous. It is not difficult to show that the operator Ski 
is entangled . A direct computation shows that 

QP(CH(Sk 3 )) = (19) 

An easy "lifting" of this construction allows to get a similar ex- 
ample for all N > 3. From the other hand , for N = 2 all rank 
nondecreasing positive operators have linear nonsingular Quantum 
perfect matchings . 

Proposition 2.7: Assuming that the Axiom of Choice and the Con- 
tinium Hypothesis hold, Ski has a Quantum semi-perfect matching 

Proof: (Sketch) Let us well order the projective unit sphere PS2 
in C 3 : S2 — {t a ; a G r) in such way that for any (3 G T the 
interval (t a : a < (3) is at most countable . Our goal is to build 
(g a ; a G F : g a / 0, g a G t„) such that if (t ai , t a2 , t a3 ) is or- 
thogonal basis then (g ai , g a2 , g a3 ) is a basis . 
As it usually happens in inductive consructions , we will induc- 
tively force an additional property : < g a ,g0 >7^ if a > (3 
and linear space L{g a ,gp) generated by (g a ,gp) is not equal to 
L(t a , tp) if < t a , tp >= 0. In this , orthogonal case , L(g a ,gp) = 
L(t a ,tp) iff g a = tp and gp — t a . Using countability assump- 
tion , it is easy to show that at each step of trasfinite induction the 
set of 'bad" candidates has measure zero , which allows always 
to choose a "good" guy p 7 without changing already constructed 
{g a ;a < 7).! 

The next Proposition shows that Shi does not have Quantum 
perfect matchings ! 

Proposition 2.8: Ski does not have Quantum perfect matchings 

Proof: Suppose that G{.) is Quantum perfect matching for Ski ■ 
We will get a contradiction by showing that then there exists a basis 
(61,62,63) such that < 61,62 >= Oand (G(6i), G(6 2 ), G(6 3 )) 
are linearly dependent . For doing that , we need to show that there 
exists an orthogonal basis (Oi, 02,03) such that O3 does not be- 
long to L(G(Oi), G(02)). Indeed , if non-zero 
de L(G(Oi),G(0 2 )) ± 
then there is no basis (G(Oi), G(02), v) with 
ved ± = L(G(Oi), G(0 2 )) , but 
(Oi , O2, d is a basis since < d, Oi . 

Take any non-zero x and an orthogonal basis {y, 2} in x 1 - such that 

G{x) = (0, 01,02) in {x, y, z} basis and 01 0, 02 =fc 0. 

LetG(y) = (61, 0, 62), G(z) = (ci,c 2 ,0) . 

Suppose that z G L(G(x), )G(y)), and y G L(G(x), )G(z)). 

Then 61 = and ci = . This contradicts 

to ((G(x), )G(y), G(z)) being a basis . Thus there exists 



an orthogonal basis (Oi, 02,03) such that O3 
does not belong to L(G(Oi), G(02)) and we got a final contradic- 
tion. ■ 

Next result shows that for weakly separable (and thus for sepa- 
rable) operators the situation is very different. 

Theorem 2.9: Suppose that T : M(N) -> M(N) is linear pos- 
itive weakly separable operator , i.e. there exists a a family of 
rank one matrices {x\y\, ...,xiyj} C M(N) such that for pos- 
itive semidefinite matrices X y the following identity holds : 

1 

Im(T(X)) = Im(^2 x iV \X yi x\) (20) 
i=i 

Then the following conditions are equivalent : 

1. T is rank non-decreasing . 

2. The rank of intersection of two geometric matroids MI(X, Y) 
is equal to N. 

3. The exists a nonsingular matrix A such that Im(AXA^) C 
Im(T(X)),X y 0. 

If , additionaly , T is completely positive then these conditions are 
equivalent to existence of nonsingular matrix A such that operator 
T'(X) = T(X) - AXA^ is completely positive . 
In this case QP{CH{T)) > N\\Det(A)\ 2 > . 

Proof: Recall Edmonds-Rado Theorem for MI(X, Y): 
Rank of MI(X, Y) is equal TV iff 

dim(L(xi;i G A) + dim(L( yj ; j G A) > N, (21) 

where Ac {1,2, 1} and A is a complement of A. 
Suppose that rank of MI(X, Y) is equal to N. Then 

RankT(X) = dim(L(a: l ; % G A)) where A =: {i : y\Xy % £ 0} 

As dim(L(2/ j; j G A) < dim(Ker(X)) = N — Rank(X) 
hence , from Edmonds-Rado Theorem we get that RankT(X) > 
N — (N — Rank(X)) = Rank(X) . 

Suppose that T is rank non-decreasing and for any Ac {1,2,...,/} 
consider an orthogonal proejctor P > on L(yj;j G A) 1 - . Then 

dim(L(x t :icA))> RankT(P) > Rank(P) = 

= N-dim(L( yj ;j CA)). 

It follows from Edmonds-Rado Theorem that rank of MI(X, Y) 
is equal to N . 

All "equivalencies" follow now directly . | 

Remark 2.10: Let us explain why Conjectures (1,2) generalize 
Hall's theorem . Consider a square weighted incidence matrix Ar 
of a bipartite graph T , i.e. Ar > if i from the first part is ad- 
jacent to j from the second part and equal to zero otherwise. Then 
Hall's theorem can be immediately reformulated as follows : A per- 
fect matching , which is just a permutation in this bipartite case , 
exists iff I Ari'|+ > |a;|+ for any vector x with nonnegative entries 
, where |a;|+ stands for a number of positive entries of a vector x . 
One also can look at Theorem(2) as a Hall's like reformulation of 
Edmonds-Rado theorem . | 



2.1 A pleminary summary 



(BUDM ) : 



So far , we got neccessary and sufficient conditions for the exis- 
tence of Quantum Perfect Matchings and presented , based on them 
, a new topological insight on the nature of the Quantum Entangle- 
ment. It is not clear to us how crucial are "logical" assumptions 
in Prop.(2.7) . Theorem(2.9) shows that in separable (even weakly 
separable) case these assumptions are not needed . The next ques- 
tion , which we study in the next sections , is about efficient , i.e. 
polynomial time , deterministic algorithms to check the existence 
of Quantum Perfect Matchings . We will describe and analyse be- 
low in the paper a "direct" deterministic polynomial time algorithm 
for weakly separable case . A complexity bound for a separable 
case is slightly better than for just weakly separable case . Our 
algorithm is an operator generalization of Sinkhorn's iterative scal- 
ing . We conjecture that without some kind of separability promise 
checking the existence of Quantum Perfect Matchings is "HARD" 
even for completely positive operators. 

3 Operator Sinkhorn's iterative scaling 

{aij : 1 < i,j < N} row 



Recall that for a square matrix A 
scaling is defined as 

R(A) = { 



J2j a a 



} 



column scaling as C( A) = { 



-} assuming that all denomina- 



tors are nonzero. 

The iterative process ...CRCR(A) is called Sinkhorn's itera- 
tive scaling (SI). There are two mainwell known properties of this 
iterative process , which we will generalize to positive Operators. 

Proposition 3.1: 

/. Suppose that A = {a itj > : 1 < i, j < N}. Then (SI) 
convergess iff A is matching, i.e., there exists a permutation 
•7T such that a ii7I (i) > (1 < i < N). 

2. If A is indecomposable, i.e., A has a doubly-stochastic pat- 
tern and is fully indecomposable in the usual sense , then 
(SI) converges exponentially fast. Also in this case there ex- 
ist unique positive diagonal matrices Di, D2, det(Z)2) = 1 
such that the matrix D^ 1 AD^ 1 is doubly stochastic. 



Definition 3.2: [Operator scaling ] Consider linear positive oper- 
ator T : M(N) — > M(N) . Define a new positive operator , 
Operator scaling , Sc 1 ,c 2 (T) as : 



S Cl ,c 2 {T){X) =: CiT(C 2 + XC 2 )C 1 t 



Assuming that both T(I) and T*(I) are nonsingular we define 
analogs of row and column scalings : 



S Ci ,c 2 {pa,b) = Ci ® C 2 {pa,b)C\ <g> C 2 f ; 

_ 1 _i 
R(pa,b) = p A 2 ® I(pa,b)p a 2 <8> I, 
_i _i 
C(pa,b) = i"® P b 2 {pa,b)I® p B 2 ■ 



(24) 



Let us introduce a class of locally scalable functionals (LSF ) de- 
fined on a set of positive linear operators , i.e. functionals satisfying 
the following identity : 

<p(Sc 1} c 2 (T)) = Det{CiCt)Det(C2Cl)ip(T) (25) 

We will call (LSF ) bounded if there exists a function / such that 

\<P( T )\ < / (M r (-0) • Ic is clear that bounded (LSF ) are natural 
"potentials" for analyzing (OSI) . Indeed , Let T n ,To = T be a 
trajectory of (OSI) , T is a positive linear operator . Then Ti (/) = I 
for odd i and T 2i (I)* = I, i > 1 . Thus if <p(.) is (LSF ) then 

9?(T i+ i) = a(i)<p(Ti),a(i) = Detail))' 1 if i is odd , 

a(i) = Det(T l (J)) _1 if i > is even. (26) 

As tr(Ti(I)) = tr(T*(I)) = N,i > , thus by the ariph- 
metic/geometric means inequality we have that | ip(T i+ \) \ > | ^(T^) | 
and if ip(.) is bounded and |<p(T)| / then DS(T„) converges to 
zero . 

To prove a generalization of Statement 1 in Prop. (3.1) we need 
to "invent" a bounded (LSF ) <p(.) such that <p(T) / iff operator 
T is matching . We call such functionals responsible for match- 
ing . It is easy to prove that QP(CH(T)) is a bounded (LSF ) . 
Thus if QP(CH(T)) =£ then DS(T„) converges to zero and , 
by Prop. (2.5) , T is rank nondecreasing . From the other hand , 
QP(CH(Sk 3 )) — and Sfc 3 is rank nondecreasing (even inde- 
composable ). This is another "strangeness" of entangled opera- 
tors , we wonder if it is possible to have "nice" , say polynomial 
with integer coefficients , responsible for matching (LSF ) ? We 
introduce below responsible for matching bounded (LSF ) and it is 
non-differentiable . 

Definition 3.4: For a positive operator T : M(N) -» M(N), we 
define its capacity as 

Cap(T) = m{{Det(X) : X y 0, Det(X) = 1} . (27) 



(22) it is easy to see that Cap(T) is (LSF ) . 



Since Cap{T) < Det{T(I)) < ( 
hence Cap(T) is bounded (LSF ) . 



tr(T(I)) jN 



R(T) = S 1 (T),C(T)=S i(T) (23) 

Operator Sinkhorn's iterative scaling (OSI) is the iterative process 
...CRCR(T) I 

Remark 3.3: Using Choi's representation of the operator T in Def- 
inition 1.5) , we can define analogs of operator scaling (which are 
nothing but so called local transformations ) and (OSI) in terms of 



Lemma 3.5: A positive operator T : M(N) — ► M(N) is positive 
rank nondecreasing iff Cap(T) > . 

Proof: Let us fix an orthonormal basis (unitary matrix) U = {u\, 
in C N and associate with positive operator T the following positive 
operator : 

Tu(X) =: ^ T(uiu\)tr(Xuiu\). (28) 

l<i<JV 

(In physics words , Tjj is a decohorence respect to the basis U , i.e. 
in this basis applying T v to matrix X is the same as applying T to 



the diagonal restriction of X. ) 

It is easy to see that a positive operator T is rank nondecreasing iff 
operators Tu are rank nondecreasing for all unitary U . 
And for fixed U all properties of Tu are defined by the following 
TV-tuple of TV x TV positive semidefinite matrices : 



A T ,u =: (T(uiu[),...,T(unv/ n 



(29) 



Importantly for us , Tu is rank nondecreasing iff the mixed dis- 
criminant M(T(uiu{), T(u N u\f)) > 0. 
Define capacity of At,u , 

Cap(A T ,u) 

ini{Det{Y Jl ^ N T{u t u\) ll ) : 7< > O^^Ji = 1}- 

It is clear from the definitions that Cap(T) is equal to infimum of 

Cap{A.T,u) over all unitary U. 

One of the main results of 1301] states that 

M(A T ,u) =: M(T(uiu\),...,T(u N ulf)) < Cap(A T ,v) < 
N N 

< -j^-M(T(uiu\), ...,T(uNU f N )). (30) 

As the mixed discriminant is a continuous (analytic ) functional and 
the group SU (TV) of unitary matrices is compact , we get the next 
inequality : 

N N 

min M(A T u) < Cap(T) < — — min M(A T u) 

uesu(N) TV! uesu(N) 

(3D 

The last inequality proves that Cap(T) > iff positive operatorT 
is rank nondecreasing. | 

So , the capacity is a bounded (LSF ) responsible for matching 
, which proves the next theorem : 

Theorem 3.6: 

1. Let T n ,To = T be a trajectory of (OSI) , T is a positive 
linear operator . Then DS(T„) converges to zero iffT is 
rank nondecreasing . 

2. Positive linear operator T is rank nondecreasing iff for all 
e > there exists e-doubly stochastic operator scaling of T 



The next theorem generalizes second part of Prop. (3.1) and is 
proved on almost the same lines as Lemmas 24,25,26,27 in [ po| ] . 

Theorem 3.7: 

1. There exist nonsingular matrices C\ , C2 such that Sc\ ,C 2 {T) 
is doubly stochastic iff the infimum in (26) is achieved . 
Moreover, ifCap(T) = Det(T(C)) where C y 0,Det(C) = 
1 

then S -11 (T) is doubly stochastic . 

Positive operator T is indecomposable iff the infimum in ( 
27) is achieved and unique . 

2. Doubly stochastic operator T is indecomposable iff 
tr(T(X)) 2 < a tr(X) 2 for some < a < 1 and all trace- 
less hermitian matrices X. 

3. If Positive operator T is indecomposable then DS(T n ) con- 
verges to zero with the exponential rate , i.e. DS(T„) < 
Ka n for some K and < a < 1 . 



4 Lower and upper bounds on Quantum Permanents 

The next proposition follows fairly directly from the second part of 
Prop.(1.9) and Cauchy-Schwarz inequality 

Proposition 4.1: Suppose that pA,B is fBUDM ). Then 



max \D(A 1 ^ (1) ,...,A N ^ (N) )\ = 
D(Ai,i, ...,Ai,n) 



(32) 



Corollary 4.2: If p A ,B is (BUDM ) then 

QP(pa,b) < NW(Ai, u ..., A ltN ) < N\Det( PA ). (33) 

Permanental part of Examplef 1. 10) shows that TV! is exact constant 
in both parts of (32) . 

The next proposition follows from the Hadamard's inequality : 
if X y is TV x TV matrix then Det(X) < Y\f =1 X(i,i). 

Proposition 4.3: IfXyO then the following inequality holds : 



Deti^TxivjXyixl) > 
1=1 

Det(X)MP (x , Y) . 



(34) 



Corollary 4.4: Suppose that separable fBUDM ) pa,b is Choi's 
representation of completely positive operator T . 
Then for all X y the next inequality holds : 

Det{T{X)) > QP{ PA , B )Det(X) (35) 

Since p A = T{I) , hence QP{p A ,b) < Det(p A ) in separable 



Call (BUDM ) p A ,B doubly stochastic if it is Choi's represen- 
tation of completely positive doubly stochastic operator T . I.e. 
(BUDM ) p A ,B is doubly stochastic iff p A = p B = I , As we 
already explained , the set of separable (BUDM ) is convex and 
closed . Thus the set of doubly stochastic separable (BUDM ) , 
DSEP(N, TV) , is a convex compact . Define 

/3(TV) = min QP(p). 

p£DSEP(N,N) 

Then it follows that (3(N) > for all integers TV . The next con- 
jecture is , in a sense , a third generation of the famous van der 
Waerden conjecture . First generation is a permanental conjecture 
proved by Falikman and Egorychev (Jl5|] , Jl4|]) in 1980 and second 
generation is Mixed discriminants conjecture po sed by R.Bapat [^] 
in 1989 and proved by the author in 1999 [|l9||. Mixed discrimi- 
nants conjecture corresponds to block-diagonal doubly stochastic 
(BUDM ) . Any good lower bound on /3(TV) will provide simi- 
larly to I pOj ] deterministic poly-time approximations for Matroidal 
permanents and new sufficient conditions for the Quantum Entan- 
glement. 



Conjecture 4.5: 



It is true for TV = 2 . | 



TV! 



(36) 



5 Polynomial time deterministic algorithm for (HMIP ) 

We introduced Hidden Matroids Intersection Problem (HMIP ) as a 
well posed computer science problem , which , seemingly , requires 
no "Quantum" background . Also , we explained that (HMIP ) 
can be formulated in terms of weakly separable (BUDM ) . Let us 
consider the following three properties of (BUDM ) pa,b ■ ( We 
will view this Pa,b as Choi's representation of completely positive 
operator T , i.e. p A ,B = CH(T) . ) 

PI Im(pA,B) contains a nonsingular matrix . 

P2 The Quantum permanent QP(pa,b) > . 

P3 Operator T is rank nondecreasing . 

We proved already that PI — > P2 — > P3 and illustrated that 
that the implication P2 — > P3 is strict . In fact the implication 
PI — > P2 is also strict. But , our Theorem (2.9), which is just an 
easy adoptation of Edmonds-Rado theorem , shows that for weakly 
separable (BUDM ) the three properties PI, P2, P3 are equivalent 
. Recall that to check PI without the weak separability promise 
is the same as to check whether given linear subspace of M (TV) 
contains a nonsingular matrix and it is very unlikely that this de- 
sision problem can be solved in Polynomial Deterministic time . 
Next , we will desribe and analyze Polynomial time deterministic 
algorithm to check whether P3 holds provided that it is promised 
that pa.b is weakly separable . 

In terms of Operator Sinkhorn's iterative scaling (OSI) we need to 
check if there exists n such that DS(T n ) < -k . If L =: min{n : 
DS(T n ) < jf} is bounded by a polynomial in TV and number of 
bits of pa,b then we have a Polynomial time Deterministic algo- 
rithm to solve (HMIP ) . Algorithms of this kind for "classical" 
matching problem appeared independently in [ p4] ] and |^] . In the 
"classical" case they are just another , conseptually simple , but 
far from optimal , poly-time algorithms to check whether a perfect 
matching exists . But for (HMIP ) , our , Operator Sinkhorn's it- 
erative scaling based approach seems to be the only possibility ? 
Assume that , without loss of generality , that all entries of pa,b 
are integer numbers and their maximum magnitude is Q. Then 
Det{pA) < (QN) N by the Hadamard's inequality. If QP(pa,b) > 
then necessary QP(pa,b) > 1 for it is an integer number. Thus 



QP{CH{T^)) = 



QP(CH(T)) 
Det(p A ) 



> (QNY 



Each nth iteration (n < L ) after the first one will multiply the 
Quantum permanent by Det(X)~ , where X >- 0, tr(X) = TV 



and tr((X - I) 2 ) > i . Using results from [E4J , Det(X)' 1 > 



(1- 



=: 5 . Putting all this together , we get the follow- 



ing upper bound on L , the number of steps in (OSI) to reach the 
"boundary" DS(T n ) < i : 



s l < QP(CH(T L )) 



(QN)- 



(37) 



It follows frm Prop.(4.2) and Cor.(4.4) that in weakly separable 

case QP(CH(T L )) < TV! 

and in separable case QP(C H(Tl)) < 1 . 

Taking logarithms we get that in weakly separable case 

L <^ 3TV(TV ln(TV) + TV(ln(TV) + ln(Q)); (38) 

and in separable case 

L <^ 3TV(TV(ln(TV) + ln(Q)). (39) 



To finish our analysis , we need to evaluate a complexity of each 
step of (OSI) . 

Recall thatT n (X) = L n {T{RlXR n ))Li , 
T n (J) = L n (T(RlR n ))Li andT*(J) = R n (T*{LtL n ))Rt . 
To evaluate DS(T n ) we need to compute tr((T£(I) —I) 2 ) for odd 
n , 

and tr((T n (I) - I) 2 ) for even n . 

Define P„ — L^Ln, Q n = R„R„ . It is easy to see that the matrix 

T n (I) is similar to P n T(Q n ) , and T* (I) is similar to Q n T* (P n ) . 

As traces of similar matrices are equal , therefore to evaluate DS(T n ) 

it is sufficient to compute matrices P n , Q n . 

But , P„+i = (T(Q„))- 1 and Q n+1 = (T^P,,))" 1 . 

And this leads to standard , rational , matrix operations with O (TV 3 ) 

per one iteration in (OSI) . 

Notice that our original definition of (OSI) requires computation of 
an operator square root . It can be replaced by the Cholesky fac- 
torization , which still requires computing scalar square roots . But 
our final algorithm is rational ! 

6 Weak Membership Problem for a convex compact set of 
normalized bipartite separable density matrices is NP- 
HARD 

One of the main research activities in Quantum Information Theory 
is a search for "operational" criterium for the separability . We will 
show in this section that , in a sense defined below , the problem is 
NP-HARD even for bipartite normalized density matrices provided 
that each part is large (each "particle" has large number of levels). 
First , we need to recall some basic notions from computational 
convex geometry. 

6.1 Algorithmic aspects of convex sets 

We will follow (l|]. 

Definition 6.1: A proper ( i.e. with nonempty interior ) convex 
set K C R" called well-bounded a-centered if there exist ratio- 
nal vector a g K and positive (rational ) numbers r, R such that 

B(a, r) C K and K C B(a, R) (here B(a, r) = {x : \\x - a\\ < 
r} and | . | is a standard eucleadian norm in R n ) . Encoding length 
of such convex set K is 

<iV>=rt+<r> + <P> + <a>, 

where <r>,<R>,<a> are the number of bits of corre- 
sponding rational numbers and rational vector . 
Following [ ]1 8[ ] we define S(K , S) as a union of all <5-balls with cen- 
ters belonging to K ; and S(K, -5) = {x G K : B(x, S) C K } . 



Definition 6.2 The Weak Membership Problem (WMEM{K, y, 8)) 
is defined as follows : 

Given a rational vector y G R n and a rational number 8 > either 

(i) assert that y G S(K, 8) , or 

(ii) assert that y g S(K, -8) . 

The Weak Validity Problem (WVAL(K, c, 7, 8)) is defined as 
follows : 

Given a rational vector y G R n , rational number 7 and a rational 
number 8 > either 

(i) assert that < c, x >=: c T x < 7 + 8 for all x G S(K, —8) , or 

(ii) assert that c T x > 7 — 8 for some x G S(K, 8) . 



In any case , L is polynomial in the dimension TV and the number 
of bits log(Q). 



Remark 6.3: Define M(K, c) =: max l6 K < c, x > . It is easy 
to see that 

M(K, c) > M(S(K, -5), c) > M(K, c) - ||c||5f ; 
M(if,c) < M(S(K,5),c) > M(K,c) + \\c\\5 



Recall that seminal Yudin - Nemirovskii theorem (tt, [|l8||) implies 
that if there exists a deterministic algorithm solving IV MEM(K, y, 5) 
in Poly(< K> + <y> + <5>) steps then there exists a de- 
terministic algorithm solving WVAL(K, c, 7, S) in Poly(< K > 
+ <c> + <(5> + <7>) steps. 

Let us denote as N SEP(M, N) a compact convex set of separable 
density matrices p A ,B : C M ® -> C* M (g> , tr{p A , B ) = 1 
, M > N . Recall that 

NSEP(M,N) = CO({a;a; t <g>j/j/ + : a; G C M ,j/ G C*^; = ||y|| = 

where CO(X) stands for a convex hull generated by a set X. 
Our goal is to prove that Weak Membership Problem for NSEP(M, N) 
is NP-HARD . As we are going to use Yudin - Nemirovskii theo- 
rem , it is sufficient to prove that WVAL(NSEP(M, N), c, 7, 8) 
is NP-HARD respect to the complexity measure (M+ < c > + < 
S > + < 7 >) and to show that < NSEP(M, N) > is polyno- 
mial in M. 

6.2 Geometry of < NSEP(M, N) > 

First , < NSEP(M, N) > can be viewed is a proper convex sub- 
set of the hyperplane in it 



The standard euclidean norm 
in it corresponds to the Frobenius norm for density matri- 

ces , i.e. \\p\\ F = tr(pp i ). The matrix j-^I G NSEP(N,N) 
and || j^zl — xx* ® J/J/ ||f < 1 for all norm one vectors x,y. 
Thus NSEP(M,N) is covered by the ball B(j±gl, 1) . Next 
we will show that B(^I, y^) C NSEP(N, N) . Recall that 
p G SEP(N, N) iff tr(CH(T)p) > for all positive operators 
T : M(i\n_- > M(N). This rather straightforward result was first 
proved in ph . Let p — {Aij : 1 < i,j < JV} be a block matrix 
as in (5). For a linear operator : M (N) — > M(N) define 
p v = {*(A iJ ):l<i,j^N}. 

The following proposition is an easy reformulation of the above 
Woronowicz's criterium . 

Proposition 6.4: p G SEP(N,N) if and only if p* y for all 

positive operators ^ such that ^(J) = I . 



Lemma 6.5: Suppose that $ : M(N) — > M(N) is linear positive 
operator and *S?(I) = 7 . 
TTjen ||*(A)|| F < v^PHf. 

Proof: For A G M(jV) denote ||A|| the operator norm induced by 
a standard euclidean norm in C (i.e. ||j4|| is the largest singular 
value of A . Recall that \\A\\ 2 < \\A\\ 2 F < N\\A\\ 2 . Let B be a 
hermitian N x N complex matrix , then \\B\\I y B y — ||B||J. 
Thus using positivity and linearity we get that \\B\\I >r ^(B) y 
-\\B\\I . We conclude that 



*(£)|| < \\ B \\ for hermitian B. 



(40) 



(The last inequality is in fact true for all matrices B ). 
Let us consider an arbitrary A G M(N) and decompose it uniquely 
as A — H\ + iff 2 where matrices H\,Hi are hermitian : 2H\ = 
A + A\ 2H 2 = -i(A - A f ). It is easy to check that 

PHI = 11^111 + 11^111. 



Therefore 

\MA)\\% = ||*(lfc)||!-+||*(ff 9 )||| < NQ^^f+lMH^f). 
By (40 ) , we get that 

||#(#i)|| 2 +||#(# 2 )|| 2 < ||Hi|| 2 +||H2|| 2 < ||£Ti \\%+\\H2 ||l = ||A|||. 
Putting all this together , we finally get that 

\\^(A)\\f < VN\\A\\ F (41) 



Theorem 6.6: Let A. be a block hermitian matrix as in (5) . If 
II A II f < V7F ^en the the block matrix I + A is separable. 

Proof: Let us consider positive linear operator ^ : M(N) — > 
Wl{N) satisfying *(/)=!. 

Then (I + A)* = I + A*. Applying inequality (41) to each 
block of A and summing all of them we get that ||A*||f < 1 . As 
the matrix A* is hermitian , we conclude that (J + A)* y 0. It 
follows from Proposition(6.4) that I + A is separable. | 
Summarizing , we get that 

B{ P'7W^ )CNSEP(N ' N)CB{ P A) 

and conclude that < NSEP(N, N) >< Poly(N). It is easy to 

get from the last inequality that < NSEP(M,N) >< Poly(m&x(N, M) 

It is left to prove that W VAL(NSEP(M, N),c,j,S) is NP-HARD 

respect to the complexity measure (M+ <c> + <8> + < 

7 >) • 



6.3 Proof of Hardness 

Let us consider the following hermitian block matrix : 
/ 

C = 



\ 






Ax . 


• Am-i 


A x 


. 





Am-i 


. 






(42) 



i.e. blocks are zero if either i 7^ 1 or j 7^ 1 and (1, 1) block 
is also zero ; Ai, Am-i are real symmetric N x TV matrices . 

Proposition 6.7: 

mn.x peNSEP(M ,N) tr{Cp) = 
m&x yeR«,\\y\\=iJ2i< l <M-i(y TA ^ 2 - 
Proof: First , by linearity and the fact that the set of extreme points 

Ext(NSEP(M,N)) = 
{xx^ ® J/j/ : x G C M ,y G C N ; \\x\\ = \\y\\ = 1} 

we get that 

max p6JVSB p(]v,jv) tr(Cp) = 
max M t® M t :!!: 60M lS eC;||x||=||H||=i tr(C(xx^ ® yj/ 1 ")). 

But tr(C(yy' ® ra')) = tr{A{y)xx' < ) , where real symmetric 
M x M matrix A(y) is defined as follows : 



A(y) 



( 

ai 
\ CbM- 



Ol 





O.M-1 




;a t =tr(A 1 yy ] ),l < i < M-l. 



Thus 

max pSJVSBP{]ViJ v) tr[Cp) = 

max|| H || =1 \ ma xA(y). 

(Above \ ma xA(y) is a maximum eigenvalue of A(y)) 

It is easy to see A(y) has only two non-zero eigenvalues (d, —d) , 

where d = Ei<»<m-i {tr(Aiyy*)f ■ 

As Ai,l < i < N — 1 are real symmetric matrices we finally get 
that 

max tr(Cp) = max \ (y T Aiy) 2 . 

p£NSEP(M,N) yeR N ,||a||=l ' 

l<i<JV-l 

I 

Proposition(6.7) and Remark(6.3) suggest that in order to prove 
NP-HARDness of 

WVAL(NSEP(M, N),c, 7, <S) respect to the complexity mea- 
sure M+ <c> + <<5> + <7>itis sufficient to prove that 
the following problem of is NP-HARD : 

Definition 6.8: (RSDF problem) Given k I x I real rational sym- 
metric matrices (Ai, 1 < i < I) and rational numbers (7,5) to 
check whether 

7 + 5 > max xeR , Alxll=1 f(x) >j-5,f(x) =: ^ {x T Aixf . 

l<i<l 

respect to the complexity measure 

(I + k + Ei< s < ; <Ai> + <6> + <i>).i 

It was shown in [^| that RSDF problem is NP-HARD provided 
k > + 1. We summarize all this in the following theorem 

Theorem 6.9: The Weak Membership Problem for NSEP(M, N) 
is NP-HARD ifN <M< WfH + 2 . 

Remark 6.10: It is easy exercise to prove that (BUDM ) pa,b 
written in block form (5) is real separable iff it is separable and 
all the blocks in (5) are real symmetric matrices . It follows that 
, with obvious modifications , Theorem 6.9 is valid for the real 
separability too . 

The construction (42) was inspired by Arkadi Nemirovski proof of 
NP-HARDness to check the positivity of a given operator [0] . | 

7 Concluding Remarks 

Many ideas of this paper were suggested by J3(i| ] . The world of 
mathematical interconnections is very unpredictable (and thus is so 
exciting) . The main technical result in a very recent breaktrough 
in Communicational Complexity [ p3| is a rediscovery of particular 
, rank one , case of a general , matrix tuples scaling , result proved 
in po| ] with much simpler proof than in [ [33| . Perhaps this our 
paper will produce something new in Quantum Communicational 
Complexity ? 

We still don't know whether there is a deterministic poly-time algo- 
rithm to check whether given completely positive operator is rank 
nondecreasing . And this question is related to lower bounds on 
Cap(T) provided that Choi's representation CH(T) is an integer 
semidefinite matrix . 

Theorem(6.9) together with other results from our paper gives a 
new , classical complexity based , insight on the nature of the Quan- 
tum Entanglement and , in a sense , closes a long line of research in 



Quantum Information Theory . Still many open questions remained 
(for the author) , for instance , is it still NP-HARD for (M, N) bi- 
partite systems wnen TV is a fixed constant ? 
We hope that the constructions introduced in this paper , espe- 
cially Quantum Permanent , will have a promising future . The 
"third generation" of van der Waerden conjecture we introduced 
above will recmire the "second generation" of Alexandrov-Fenchel 
inequalities ||lj. We think , that in general , mixed discriminants 
and mixed volumes should be studied (used ) more enthusiastically 
in the Quantum context . After all , they are noncommutative gen- 
eralizations of the permanent .... 

Most of all , we hope that a reader will be able to "factor" our lousy 
english and to see the subject . 

It is my great pleasure to thank myLANL colleagues Manny Knill 
and Howard Barnum . 

Finally , I would like to thank Arkadi Nemirovski for many enlight- 
ening discussions . 
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